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ABSTRACT 



A small subset of thermonuclear X-ray bursts on neutron stars exhibit such a strong photospheric expansion that for a few seconds the 
photosphere is located at a radius r p h i 10 3 km. Such 'superexpansions' imply a large and rapid energy release, a feature characteristic 
of pure He burst models. Previous calculations have shown that during a pure He burst, the freshly synthesized heavy-element ashes 
of burning can be ejected in a strong radiative wind and produce significant spectral absorption features. We search the burst data 
catalogs and literature and find 32 superexpansion bursts, 24 of which were detected with BeppoSAX and three with RXTE at high 
time resolution. We find that these bursts exhibit the following interesting features: (1) At least 31 are from (candidate) ultracompact 
X-ray binaries in which the neutron star accretes hydrogen-deficient fuel, suggesting that these bursts indeed ignite in a helium-rich 
layer. (2) In two of the RXTE bursts we detect strong absorption edges during the expansion phase. The edge energies and depths are 
consistent with the H-like or He-like edge of iron-peak elements with abundances ^100 times solar, suggesting that we are seeing 
the exposed ashes of nuclear burning. (3) The superexpansion phase is always followed by a moderate expansion phase during which 
r ph ~ 30 km and the luminosity is near the Eddington limit. (4) The decay time of the bursts, T decay , ranges from short (« 10 s) 
to intermediate 10 3 s). However, despite the large range of Tdecay* me duration of the superexpansion is always a few seconds, 
independent of T decay . By contrast, the duration of the moderate expansion is always of order T it . c . dy . (5) The photospheric radii r ph 
during the moderate expansion phase are much smaller than steady state wind models predict. We show that this may be further 
indication that the wind contains highly non-solar abundances of heavy elements. 

Key words. X-rays: binaries - X-rays: bursts - Nuclear reactions, nucleosynthesis, abundances - stars: neutron - radiative transfer - 
X-rays: individual (4U 1722-30, 4U 0614+09, 4U 1820-30) 



1. Introduction 

X-ray bursts are thermonuclear shell flashes that result from the 
unstable burning of fuel accreted on the surface o f a neutron star 
(NS) f rom a Roche-lobe filling c ompanion star (|Grindlav et alj 
197 6L IWoosle v & Taaml 1 19761 iMaraschi & Cavalierel 11977: 
lLamb & Lambl 119781: for reviews^ see iLewin et al.l 11995 



Strohmayer & Bildsten 2006). Approximately 20% of all bursts 
exhibit s ome degree of photosp heric radius expansion (PRE; 
see, e.g., lGallowavetaiT l2008a) due to a flux exceeding the 
Eddington limit, res ulting in a radiat i vely driven expansion 
of the photosphere dLewin et all Il984t lEbisuzaki etail 119831: 
, iTawara et al.lll984allbl) . During the expansion, the emission area 
increases while the bolometric luminosity remains nearly con- 
stant (near the Eddington limit). As a result, the photosphere 
cools and the detected X-ray flux decreases. Eventually the bolo- 
metric luminosity decreases, the photosphere returns to the NS 
surface, and the X-ray flux rises again, signifying the start of the 
main burst. 

In a small fraction of bursts (we estimate no more than a few 
tenths of a percent), the PRE is so extreme that for a few sec- 
onds the detected X-ray flux decreases to less than 1 % of its peak 
value. Such extreme PREs cre ate the illusion of a pre cursor to the 
main burst, as first reported bv lHoffman et al.l(ll978l) . Given that 
the peak temperature of a burst is typically 2.5 keV, close to the 
soft end of a typical X-ray device's bandpass, such a dip implies 
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that the effective temperature decreases to less than 0.2 keV (as- 
suming pure black body radiation at a constant bolometric flux 
and a hydrogen column density of 10 22 H-atoms cm 2 , typical 
for an X-ray burst source). This corresponds to a radius expan- 
sion factor of 100 or more, implying that the photosphere moves 
out to radii r p h ^ 10 3 km. We call such an extreme PRE a super- 
expansion to distinguish it from the much more common mod- 
erate PRE with an expansion factor of a few. 

In X-ray burst models, bursts that ignite in helium-rich lay- 
ers have high peak luminosities that often exceed the Eddington 
limit and drive strong r adiative winds jHanawa & Sugimoto 
1982; iPaczvnskil 119831: iFuiimoto et all 119871: IWooslev et al.l 
2004). Superexpansion bursts are the most luminous of all X- 
ray bursts and may be the result of nearly pure helium ignition; 
as such they serve as a useful test bed for helium burst models. In 
addition, calculations show that heavy-element ashes of nuclear 
burning can be ejected in the winds of PR E bursts and pro duce 
strong spectral absorption features (Weinberg et al. 2006b). The 
brightness and powerful, long-lasting winds of superexpansion 
bursts makes them ideal candidates for detecting spectral evi- 
dence of heavy-element ash ejection. 

Helium bursts can be of either short or long duration de- 
pending on the properties of the accretion (Fuii moto et al.l l 1981; 
Bilds tenl[T998l) . If the NS accretes a H/He mixture, the bursts 
tend to be short (seconds). If, however, the system is an ul- 
tracompact X-ray binary (UCXB) and the NS accretes hydro- 
gen deficient material, the bursts tend to be much longer (sev- 
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eral tens of minutes , i.e., the 'intermedia te-duration bursts'; 
lin 'tZand et al 1 l2005t ICumming et aD l2006). This is because in 
UCXBs, unlike in mixed H/He accretors, there is no stable CNO 
burning to keep the freshly accreted envelope hot. As a result, 
the helium fuel ignites at a greater depth, where the thermal time 
is longer. As we will show, the superexpansion bursts that we 
find are all associated with (candidate) UCXBs. 

The volume of data on superexpansion has increased sub- 
stantially since the last detailed publication on extreme PREs 
(Ivan Paradiis et al.l T1990). At that time only 4 superexpansion 
bursts were known. It therefore makes sense to revisit the sub- 
ject. In §[2] we report on an additional 28 superexpansion bursts, 
most of them found in archival BeppoSAX and RXTE data. In 
§ [3] we examine the spectroscopic features of the best-covered 
superexpansion bursts. We find evidence for a variable absorp- 
tion edge that may be due to exposed heavy-element ashes of 
nuclear burning. In §[4] we describe the spectroscopic features 
of the other superexpansion bursts. We find that in all 32 bursts 
the superexpansion phase is followed by a moderate expansion 
phase. In §[5] we show that the duration of the two phases are 
uncorrelated. We place these results in the context of models of 
helium-rich X-ray bursts and radiation-driven NS winds in § [6] 
and conclude in § Q 

2. Search for superexpansion bursts 

We searched for bursts with precursors in the literature 
(Gall oway et alj|2008at iKuulkers et a!1l2003l and other relevant 
literature; see Table Q]) and data of the Proportional Counter 
Array (PC A; operational since January 1996) on RXTE up 
to and including 2008 and in all data of the two Wide Field 
Cameras (WFCs; operational between June 1996 and April 
2002) on BeppoSAX. Both instrument packages consist of pro- 
portional counters with similar efficiency curves and spectral 
resolution (20% full-width at h alf maximum), and bandpasses 
that start at 2 keV. The PCA dJahoda et all 120061) consists of 
5 proportional counter units (PCUs) that each have a propane- 
filled and a xenon-filled gas layer. The xenon layer is the main 
detector and is read out at high time and energy resolution 
between 2 and 60 keV. The propane layer is used for anti- 
coincidence. It has only 0.125 s time resolution and no energy 
resolution. However, since it has additional sensitivity between 
1 .5 and 2 ke V, it is useful for studying the soft phases of X-ray 
bursts. The net collecting area of the xenon layers combined is 
about 8000 cm 2 , the field of view 2° (full-width to zero response) 
and there is no angular resoluti on. The time resol ution can be set 
to 1 [is at best. The two WFCs dJager et al.ll 19971) each had a net 
collecting area of 140 cm 2 , a field of view of 40°, an angular res- 
olution of 5' in a bandpass of 2 to 28 keV. The time resolution 
was 0.5 ms. The observation program of the PCA is in large part 
dedicated to burstin g NSs, yielding large exposure times (of or- 
der megaseconds: lGallowav et aLll2008al) despite the small field 
of view. Thanks to the wide field of view of the WFCs, long ex- 
posure times were acquired on the majority of bursting NSs, re- 
sulting in an unprecedented ability to detect rare and unexpected 
events such as superexpansion bursts. 

Our search through ss 3400 bursts and the literature turned up 
32 superexpansion bursts from 8 different systems, see Table Q] 
All the bursts that have been unambiguously identified with an 
ob ject (7 out of 8 systems ) are from (candidate) UCXBs as listed 
by lin 't Zand etail (120071) . The UCXB candidacy of the only 
transie nt, XB 1715-321, has been put forward in iJonker et alJ 
J2007h . The one burst not unambiguously identified with an ob- 
ject is from 4U 1708-23. Of the 32 bursts listed, 22 are from 
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Fig. 1. Time profiles of the 2-30 keV photon flux of 24 X-ray 
bursts detected from 4U 1722-30 with BeppoSAX WFC (this 
includes 21 bursts with precursor, see Table [1] and three without 
at MJD 50521.804, 50527.181 and 50532.41 1). The photon flux 
is scaled between -1 and 4 phot s^'cirT 2 . The horizontal dotted 
lines indicate a flux level of phot s^'crrr 2 . The MJD indicated 
in each panel refers to t = s. The profiles have been aligned to 
the start time of the precursor. 



4U 1722-30 (in the globular cluster Terzan 2), 21 of which were 
detected with BeppoSAX- WFC (see Fig. IT1 and iKuulkers et al.l 
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Table 1. Measurements of 32 bursts with superexpansion. 





Duration (s)° 




Instr./MJD 




Super 


Moderate 


^"decay 


Ref.* 




Precursor 


expansion 


expansion 






phase 


phase t x 


phase f me 




4U 0614+09 


RXTE/5 1944.903 


0.05 


1.2 


100 


40(2) 


1 


A 1246-588 


WFC/50286.290 


3.0 


1.5 


54 


38(5) 


2 


WFC/5 1539.874 


1.5 


6.0 


25 


19(1) 


2 


4U 1708-23 (probably) 


SAS-C/43 18 1.834 


4.2 


6.0 


304 


300(50) 


3,4 


XB 1715-321 


SAS-C/42957.620 


2.4 


1.5 


36 


30 A 


3 


Hakucho/45 170.231 


3 


4.0 


105 


85(5) 


5 



4U 1722-30 



WFC/503 18.279 


4.0 


2.0 


16 


18(3) 


6 


WFC/50330.196 


2.0 


5.0 


20 


16(6) 


6 


WFC/50348.938 


3.0 


3.0 


15 


16(5) 


6 


WFC/50368.307 


2.0 


5.0 


14 


19(3) 


6 


WFC/50526.311 


1.0 


4.0 


15 


21(6) 


6 


WFC/50536.895 


2.5 


3.5 


11 


9(3) 


6 


WFC/50538.439 


3.0 


1.0 


15 


19(6) 


6 


WFC/50553.130 


2.5 


1.5 


17 


12(2) 


6 


WFC/50892.706 


3.0 


2.5 


18 


15(2) 


6 


WFC/50904.813 


2.5 


4.0 


21 


19(3) 


6 


WFC/5 1057.579 


1.5 


3.0 


22 


19(4) 


6 


WFC/5 123 1.379 


2.0 


5.0 


17 


26(11) 


6 


WFC/5 1270.560 


2.0 


5.5 


25 


14(2) 


6 


WFC/5 1278.690 


2.5 


2.5 


15 


13(1) 


6 


WFC/5 1422.838 


2.0 


6.5 


22 


28(5) 


6 


WFC/5 143 1.282 


4.0 


5.0 


27 


26(5) 


6 


WFC/5 1453.377 


1.5 


5.5 


20 


24(6) 


6 


WFC/51461.331 


4.0 


3.5 


19 


23(4) 


6 


WFC/51610.000 


3.0 


3.5 


20 


17(3) 


6 


WFC/5 1639.966 


1.5 


5.0 


32 


23(5) 


6 


WFC/5 1956.091 


2.0 


3.5 


41 


29(10) 


6 


RXTE/50395.292 


3.6 


1.6 


23 


30.2(0.1) 


7 


SLX 1735-269 


I'GRAL/52897.733 


2.0 


7.0 


482 


600(100) 


8 


4V 1820-30 


RXTE/5 1430.074 x 


15.0 


2.3 


1400 


2500 


9 


M15 X-2 


Ginga/47454.730 


1.5 


5.5 


88 


60 


10 


WFC/5 187 1.593 


1.5 


7.5 


169 


155(11) 





'See Fig. [4] for a plot of / sc and f m 
theses represen t lcr uncertain ties; * 
lin 't Zand et ail d2008h . 3 - iHoffman et al.1 
d!984l). 



Tdecay ^Numbe rs in paren- 
Kuulke rs et all d2009l) . 2 
dl978h. 4 



Tawaraetal. (1984b), 6 - Kuulkers et al. 



Lewin et al. 
d2003h. 7 - 



Molk ov et al1d2000h. 8 -iMolkov et alj d2005h . 9 - IStrohmaver & BrowrJ 
d2002h . 10 -Ivan Parad iis et al.1 ((1990); A This number is rather uncertain 
due to incomplete coverage of the burst; x superburst 



l200l . Almost all bursts from 4U 1722-30 have a precursor; for 
a burster with precursors, it has relatively frequent bursts (once 
every few da ys). The superexpansion bur st from 4U 1820-30 is 
a superburst (IStrohmaver & Brownll2002l hereafter SB02). 

3. Spectral analysis of the 3 PCA bursts 

The data for the three PCA superexpansion bursts are of higher 
quality than the WFC burst data due to the PCA's large pho- 
ton collecting area. The quality is especially good because at 
least three PCUs were employed during these observations (this 
is generally true of only « 5% of the PCA's observation time). 



Nonetheless, the spectroscopic resolution of the PCA is modest 
(17% FWHM at 6 keV) and only relatively broad spectral fea- 
tures can be resolved. 

High-time-resolution spectrally-resolved data are available 
for the first 30 s of the burst from 4U 0614+091. Such data are 
available for much of the 4U 1722-30 burst, except between 3 
and 8 s after burst onset. The superburst from 4U 1820-30 had 
good spectral coverage throughout but only at a modest time res- 
olution of 16 s (the 'Standard-2' data collecting mode). We de- 
scribe the spectral data and analysis method in more detail in 
Appendix lAl 

Typically burst spectra are well fit by an absorbed 
black body model (ISwank et al.1 119771; IHoffman et all 119771; 
IStrohmaver & Bildstenl 20061) . Deviations from this model have 
been reported for the s uperexpansion bursts from M15 X-2 
(Ivan Paradiis et al.l 1 1990b and the PCA superexpansion super- 
burst from 4U 1820-30 (SB02). In the latter case, SB02 detected 
a broad emission line between 5.8 and 6.4 keV and an absorption 
edge at 8 - 9 keV at times > 100 s after burst onset (at earlier 
times the blackbody temperature was evolving too quickly for 
the Standard-2 accumulations to resolve). To model these fea- 
tures, SB02 included a Gaussian emission line and an absorp- 
tion edge in their spectral model. With these components, they 
found acceptable fits with reduced x 1 ~ 1 \ both the line and 
edge were strongly required to adequately model the data. SB02 
further noted that features of the Gaussian line — its finite width 
and an energy centroid that is often significantly less than 6.4 
keV — suggested that the line was produced by reflection from a 
relativistic accretion disk. Indeed, excellent fits to the data were 
obtained by modeling the line as an Fe fluorescence feature from 
an accretion disk and accoun ting for smearing of the absorption 
edge due to disk motions dBallantvne & St rohmaver 120041 and 
lBallantvnell2004l) . 

3.1. Absorption edges 

As shown in Appendix iBl (see top row of Fig. IB. Il l, the bursts 
from 4U 0614+091 and 4U 1722-30 also show significant devi- 
ation from an absorbed black body. The deviations are similar to 
those seen in the 4U 1820-30 superburst: a broad emission line 
at 6 to 7 keV and an edge at « 10 ke V that varies with time. 

M otivated by the success of Ballan tvne & Strohmaverl 
(2004) at fitting the deviations in 4U 1820-30 with a model that 
includes a reflection component, we attempted to fit the spec- 
tra of the 4U 0614+091 and 4U 1722-30 bursts using the same 
black body plus reflection model as Ballantyne & Strohmayer 
(provided to us by D. Ballantyne, see Appendix IBl. Although, 
as shown in Fig. IB. II the model fits some of the spectra, it fails 
to fit all the spectra, particularly spectra just after the end of the 
superexpansion phase. Examples of when the black body plus 
reflection model fail are shown in the top row of Fig. [2] These 
include two spectra of 4U 0614+091 within a few seconds after 
the superexpansion and one spectrum of 4U 1722-30 during the 
first 7 s after the superexpansion (as indicated by the grey bar in 
Fig.fSJl. In these three examples, the edges are deeper or at lower 
energies than the reflection model can accommodate. 

It is worth noting that during the superexpansion, the flux 
of both bursts drops below the pre-burst level (see top panel 
of Fig. |3). This may be because the expanding shell and super- 
Eddington flux disrupt the disk and temporarily stop the accre- 
tion. Alternatively, the disk is not disrupted and we just do not 
see the accretion flux because the expanding shell surrounds and 
obscures the inner disk. Whatever the cause, it is perhaps not 
surprising that the disk reflection model fails at times immedi- 
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4U 061 4+09 at 1 .7 s / Black body + reflection fit 4U 061 4+09 at 2.5 s / Black body + reflection fit 



4U 1 722-30 / Black body + reflection fit 




4U 061 4+09 at 1 .7 s / Fit with edge included 



4U 0614+09 at 2.5 s / Fit with edge included 



4U 1722-30 / Fit with edge included 
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Fig. 2. Spectral fits to the PCA bursts from 4U 0614+09 and 4U 1722-30 during time intervals when the absorbed black body plus 
reflection model fails to fit the data (see Fig. IB. II for example time intervals when the black body plus reflection model provides a 
good fit to these bursts). Left plots: The 4U 0614+09 spectrum between 1.65 and 1.78 s after burst onset (this is the second data 
point after the superexpansion phase in Fig. [3]). Middle plots: The 4U 0614+09 spectrum between 2.53 and 2.65 s after burst onset 
(ninth data point). Right plots: The 4U 1722-30 spectrum for the 7s long interval indicated by the grey bar in Fig. [3] The upper 
panels show the observed photon flux (crosses) and the fitted black body plus reflection model (histogram). The middle panels show 
the fractional deviation of the data points from the black body plus reflection model. The lower panels show the deviations after an 
absorption edge is included in the model. 



ately following the superexpansion. Furthermore, even when the 
reflection model does fit the data, it often yields unreasonable 
results; for instance, it requires that the majority of the flux be in 
the reflection component rather than the direct black body com- 
ponent. This results in unreasonable values for the black body 
luminosity and radius. 

The failure of the reflection model to fit some of the spec- 
tral data from 4U 0614+091 and 4U 1722-30 prompted us 
to add another component to the spectral model: an absorp- 
tion edge stemming from the pr imary (non-reflected) s pheri- 
cally symmetric burst emission. IWeinberg et all (|2006b) have 
suggested that such a signature should appear in the spectra 
of radius expansion bursts due to the presence of heavy ashes 
of nuclear burning in the photosphere (see § I6.21 i. To model 
such an edge, we multiply the spectral energy distribution by 
M(E) = exp[-T opt (£/£ e dge)~ 3 ] for photon energies E > E edge 
and M(E) — 1 otherwise^]. The fit parameters are the optical 
depth of the line r opt and the edge energy £edge- Our full spectral 
model thus accounts for three features: an absorbed black body, 
the reflection of the burst flux from an inner accretion disk, and 
an absorption edge formed in the photosphere. 

We find that including such an edge significantly improves 
the fitfl as illustrated in the lower row of plots in Fig. [2] The evo- 
lution of the absorption edge parameters T opt and £edge is shown 



1 This is multiplicative model edge in XSPEC. 

2 For some spectra the edge can be satisfactorily fit by allowing for 
unrealistic values of Nu, the black body temperature and the accretion 
flux. However, for the 4U 0614+091 burst, the edge immediately after 
the superexpansion has such a low energy (see below) that a good fit 



in Fig.[3]for the three PCA bursts. The absorption edge from 4U 
0614+091 exhibits a peculiar increase in £edge from 4.63 ± 0.05 
keV just after the superexpansion to 8.5 ±0.1 keV less than 2 s 
later (see also left and middle panels of Fig. |2). For the next 2 
s, ^edge remains nearly constant at 8.5 keV, after which no sig- 
nificant edge is detected. While Sedge is increasing, T opt is large 
(between 2 and 3), after which it becomes substantially less than 
1. 

The fit to the first spectral data point from 4U 1722-30 (8 s 
after burst onset) yields an absorption edge with £ e dge — 1 1 keV 
and r opt =a 1 . The parameters remain near these values until the 
temperature peaks at 23 s after onset, after which the parame- 
ters quickly change to £ e dge -6 — 7 keV and T opt - 0.5. Since 
there is no spectral data immediately following the superexpan- 
sion (from 3 to 8 s), we do not know if the increase in Ziedge 
seen in the spectra from 4U 0614+091 occurs in the case of 4U 
1722-30. 

The absorption edge from 4U 1820-30 is shallow, but still 
detected with significance thanks to the high quality signal pro- 
vided by the long duration of the burst. This is particularly clear 
when co-adding all the data accumulated over a long time span 
(e.g., all data between 5000 and 10000 s). 

Although we obtain reasonable fits to the data once we in- 
clude an edge in the model, the limited spectral resolution of 
the PCA precludes an unambiguous identification of these fea- 
tures as absorption edges. Higher spectral resolution observa- 



cannot be obtained even if one allows for unrealistic values of these 
standard spectral parameters. 



4U 0614+091 
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4U 1722-30 
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Fig. 3. Time pr ofiles of various spec tral parameters for 3 superexpans ion bursts detected with RXTE-PCA from 4U 0614+09 
(left; see also IKuulkers et alj 120091). 4U 1722-30 (center ; see also iMolkov et all 120001) and 4U 1820-30 (right; see also 
Strohmaver & Brown 2002; Ballantvne & Strohmaver 2004). Time is in seconds since burst onset (see Table IA. lb . Top panels: 
total PC A countrate for the xenon layers (solid line, in units of 10 4 s _1 ) and propane layer (dot-dashed grey line, in units of 
1,250 s _1 ), after subtracting the pre-burst levels. The horizontal dashed lines mark a count rate of zero. Second panels: bolometric 
flux of the fitted black body if no reflection component is included, corrected for detector dead time and for absorption. The errors 
were not determined, but they become large when kT < 1 keV. Third panels: black bo dy temperature in ke V. Fourth panels: black 
body radius assumi ng isotropic emission. Assu med distances: 3.2 kp c (4U 06 14+09: IKuulkers et al.ll2009l) . 9.5 kpc (4U 1722-30; 
Kuul kers et al.ll2003h and 7.6 kpc (4U 1820-30; Kuu lkers et al.ll2.Q03h : Fifth and sixth panels: absorption edge energy Eedge in keV 
and optical depth T opt . Bottom panels: reduced/^ 2 value for each fit to the black bo dy + reflection + edge model. W e left out the data 
for 4U 1820-30 between 2000 and 2500 s, because these are not understood (see lBallantvne & Strohmaverl2004l) . The grey vertical 
bars for 4U 0614+09 and 4U 1722-30 delimit the time intervals for the spectra shown in Fig. [2] A^h was fixed for 4U 0614+09 and 
4U 1820-30 and free for 4U 1722-30 because a fixed value did not yield good fits. 



tions of superexpansion bursts, for instance with XMM-Newton 
or Chandra, would be extremely useful. 

3.2. Other spectral features of the 3 PCA bursts 

The evolution of some of the other spectral parameters is shown 
in Fig. [3] We do not show the results for the reflection parame- 
ters as these are not of immediate interest. The top panel shows 
the light curve for both the xenon and propane detector layers. 
The most conspicuous feature is the significant dip in the light 
curve during the superexpansion phase; the flux drops below the 
pre-burst level signifying that the large burst flux temporarily 
disrupts or obscures the accretion on to the NS. The propane 
layer is most sensitive at 1.5 to 2 keV and the signal from this 



layer is strongest when the black body temperature kT us 1 keV 
immediately following the superexpansion. The light curve from 
4U 1722-30 has a 2 s long spike near the end of the superexpan- 
sion phase that is especially pronounced in the propane layer. It 
may be due to a brief contraction and re-expansion of the photo- 
sphere, although interpreting the spike is difficult because there 
is no spectral data available during this time interval. There are 
indications of similar spikes in the 4U 1722-30 bursts detected 
with the WFC (Fig. [TJ. A similar feature has been reported in a 
burst from M 1 5 X-2 Ivan Paradijs et al.fl990l) . 

The 3 PCA bursts show a phase of moderate expansion after 
the superexpansion. During the moderate expansion phase, the 
black body temperature increases and the radius decreases. We 
define the end of the moderate expansion phase to be the time 
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at which the temperature peaks, which is attribute d to the pho- 
tosph e re 'touching-down' on the NS surface (e.g.. [ Lewin et al.l 
119931: iKuulkers et ail 120031; iGallowav et al . 2008a,tJ). This is 
nicely illustrated in Fig. [3] in the bursts from 4U 1722-30 and 
4U 1820-30, where touch down occurs 23 and 1400 s after burst 
onset, respectively. The moment of touch down in 4U 0614+09 
was not covered by the P CA, but was observed with FREGATE 
on HETE-II to be at 90 s (IKuulkers et alj[2009b . 

In the bursts from 4U 0614+09 and 4U 1722-30, we find 
that at the start of the superexpansion phase, when the X-ray 
flux begins to drop, the photosphere is apparently moving out at 
velocities u p h ~ 3 x 10 4 kms _1 and u p h as 10 3 kms _1 , respec- 
tively (left and middle panels of Fig. [5J this estimate is simply 
the change in apparent radius divided by the time over which that 
change occurred). It is difficult to determine whether these mea- 
surements reflect the true velocities of the photosp here due to un- 
certainties, e.g., in the bla ck body color corr ection (lLondon et al.l 
Il986t iPavlov et al.lll99lT) . Calculations by iMadei et al.l d2004) 
show that the observed color temperature of a black body fit 
may differ by a factor of 1.3 to 1.4 from the effective temper- 
ature for a 1 keV black body and about 1.6 for a 3 keV black 
body. Consequently, inferred radii would need to be corrected 
by a factor of 1 .7 to 2.6. In addition, the PCA calibrated data do 
not cover photon energies below 3 keV and, thus, the peak of the 
thermal spectrum when the expansion is large. We nonetheless 
consider the expansion velocities to be reliable to within a factor 
of a few. 

The radii during the moderate expansion phase of the 4U 
1820-30 superburst are unreasonably small (i.e., < 10 km). 
It cannot be attributed to an underestimated val ue of the dis- 
tance as this is well measured both optically (Heasle v et al.1 
2000) and from comparing the burst peak flux to the Eddington 
limit of a hydrogen-poor p hotosphere (e.g., IZdziarski et al.l 
I2007t IGallowav et al.ll2008al) . A color correction would, how- 
ever, increase the measured radii to acceptable values (see also 
IStrohmaver & Brown 2002). The superburst also behaves some- 
what unusually between 2000 and 2500 s after burst onset. The 
behavior is unexplained and the data in this time frame is ex- 
cluded in both our analysis and that of Ball a nt vne & Strohmaverl 
(200|). 



4. Spectroscopic features of the other bursts 

When there is sufficient sensitivity, other superexpansion bursts 
show some or all of the features seen in the 3 PCA bursts. The 
edge is difficult to detect in bursts measured with instruments 
other than the PCA, although the burst from Ml 5 X-2 measured 
with Ginga showed indicatio ns of edges in the spectrum (see Fig. 
6 in lvan Paradiis et aLlll990l) . The one feature that is always seen 
is a phase of moderate expansion following the superexpansion. 
The spike in low-energy flux a few seconds after superexpansion 
is not always seen (e.g., it is not seen in the PCA bursts from 4U 
0614+09 and4U 1820-30). 

The burst from 4U 0614+09 has a very fast rise and the 
precursor is merely 0.05 s long (Fig. [3J. For such short dura- 
tions, smaller detectors might fail to catch the precursor and in- 
stead first see the burst during the end of the superexpansion 
phase, when the photospheric radius is decreasing. An example 
of this might be the long burs t detected from 2S 0918-549 with 
the BeppoSAX WFCs din 't Zand et al.ll2005l) . This burst lacks a 
precursor and yet spectroscopic analysis of the rising phase re- 
veals a black body radius that is rapidly decreasing, as if the burst 
was preceded by a superexpansion phase. The very long burst 
from SLX 1737-282 (T de ca y = 659+32 s, see also lin 'tZandetalJ 
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Fig. 4. Duration of radius expansion phases versus decay time 
Tdecay for each superexpansion burst. For data, see Table [TJ 
Crosses mark the duration of a burst's superexpansion phase (f se ) 
and asterisks mark the duration of its moderate expansion phase 
(?me)- The dashed line marks the relationships f se , f me = Tdecay 



I2OO2I: iFalanga et alj|2008h also lacks a precursor despite having 
a radius that decreases during the rise, although the decrease is 
not as extreme as in 2S 0918-549. 

There are a few bursts that also show two expansion phases, 
but whose expansion factor in the first phase is too small (be- 
tween 5 and 10) to qualify as a superexpansion. Their X-ray sig- 
nal is not completely lost, but their light curve still shows a deep 
dip. This is true of all o rdinary X-ray bursts d etected from 4U 
1820-30 with the PCA (IGallowav et alj|2008al) . a second PCA 
burst from 4U 1722-30 detected at M JD 54526.679 (Galloway, 
priv. comm.), a burst from M15 X-2 dSmaldl200lh . and a few 
bursts from the non-UCXBs 4U 1705-44, KS 1731-260, SAX 
J1747.0-2853, and SAX J1808.4-3658. 



5. Measurement of expansion durations 

We define the duration of the superexpansion phase, f se , as the 
time between the end of the precursor and the start of the main 
burst. In practice (allowing for statistical uncertainties), this is 
the time over which the flux is less than about 10% of its peak 
value. We define the duration of the moderate expansion phase, 
f me , as the time between the end of the superexpansion and the 
time of peak temperature or spectral hardness, which is usually 
considered to be the moment when the expanded photosphere 
'touches down' on the NS (see §[3]). The decay time Tdecay is de- 
termined by fitting an exponential function to the full-bandpass 
light curve after touch down and measuring the e-folding decay 
time. Since all the detectors are of a similar type (i.e., a xenon- 
filled proportional counter), the systematic errors introduced by 
using observed photon rates instead of derived calibrated energy 
fluxes are negligible compared to statistical errors. 

Table[TJgives the duration measurements for the 32 superex- 
pansion bursts. The durations f se and f me are plotted as a function 
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of Tdecay in Fig. [4] The diagram reveals two important features: 
first, f me is proportional to Tdecay (they are approximately equal). 
This is not surprising; the longer the burst, the longer the flux re- 
mains near its peak value. Since the peaks are super-Eddington 
in these bursts, the duration of the super-Eddington phase should 
be proportional to the decay time (see AppendixP. Second, f se 
is independent of Tdecay; in all bursts the superexpansion lasts a 
few seconds. We discuss possible explanations for the lack of 
correlation between f se and Tdecay in § 16.31 

6. Discussion 

6.1. Association between superexpansion bursts and UCXBs 

All the superexpansion bursts that we find, in cluding those in 
the lit erature, are from (candidate) UCXBs (see lin 't Zand et aTl 
l2007l) . A cha racteristic of UC XBs is that they are hydrogen de- 
ficient (e.g.. lNelerriansl F2008). As a result, X-ray bursts from 
UCXBs are most likely fueled entirely by helium and/or heav- 
ier element^. Furthermore, since there is no stable CNO burning 
in-between bursts, the freshly accreted NS envelope in an UCXB 
is colder than that of a hydrogen-rich system accreting at the 
same rate. Bursts from UCXBs therefore tend to occur at greater 
depth, where the thermal time is longer. This gives rise to so- 
called 'intermediat e-duration bursts' w i th durations of severa l 
tens of minutes (in 't Zand et a I] l2005t [Cumming et al. 2006). 
Because of the large and rapid energy release, the peak lumi- 
nosity of helium-rich bursts (of both the short- and intermediate- 
duration variety) often exceeds the Eddingto n limit, resulting 
in phases of photospheric radius expansion {Paczvhski 1983; 
IWeinberg et al.ll2006bHCumming et al1l2006l) . 

6.2. Evidence for heavy-element ashes in photosphere 

The convective region that forms during an X-ray burst is well- 
mixed with the freshly synthesized ashes of nuclear burning. 
IWeinberg et al.l d2006 b) solve for the evolution of the vertical ex- 
tent of the convective region in PRE bursts with strong radiative 
winds. They find that the convective region extends outwards to 
sufficiently shallow depths that the base of the wind lies well 
within the ashes of burning. As a result, some ashes are ejected 
in the wind and exposed at the photosphere. The column den- 
sity of ejected ashes, which consist primarily of heavy elements 
(near the Fe peak; see below), is expected to be especially large 
for the energetic, helium-rich bursts that exhibit superexpansion. 
The spectral signature of the ashes is therefore expected to be 
detect able at e ven modest spectr al r esolution (IWeinberg et al.l 
2006b} see also lvan Paradiislll982l and lFoster et al.lll987l) 

We now discuss two features in the data that suggest the 
presence of heavy-element ashes in the photosphere: absorp- 
tion edges at the end of the superexpansion phase (§ 16.2. U 
and smaller-than-expected radii during the moderate expansion 
phase (§ I6.2.2| |. 

6.2.1. Absorption edges 

Which elements might be responsible for the observed absorp- 
tion edges? We first consider the early-time edge in the spec- 
trum of 4U 0614+09, which increases from an initial energy of 

3 It is unclear whether free-falling a particles, in the absence of pro- 
tons, will spall heavier accreted nuclei already present in the NS atmo- 
spher e to such a degree th at significant amounts of protons are produced 
(cf jBildsten et all 1991120031) . 



£edge = 4.6 keV to 8.5 keV within 2 seconds. Since the pho- 
tosphere is at large radii (r p h ^ 30 km) during this time, the 
gravitational redshift should be negligible (z ^ 0.1). Although 
r p h is decreasing with time, suggesting that material is moving 
away from the observer and back on to the NS, the rate of de- 
crease appears much too slow (~ 100 km s ) to produce any 
appreciable doppler redshiftQ This suggests that the variations 
in Sedge are not due to variations in redshift but rather are intrin- 
sic to the source. One possibility is that the shift in £edge is due 
to variations in which elements and ionization states dominate 
the absorption as the photosphere contracts and T e ff increases. 
If correct, bound-bound transitions should also be a source of 
absorption. The limited spectral resolution of the PCA makes it 
difficult to detect narrow X-ray absorption lines, however (such 
lines have never been confirmed in PCA spectra as far as we 
know). 

We now consider the late-time 8.5 keV edge in 4U 0614+09. 
The measured photospheric radii at this time (r p h ~ 10 km) im- 
ply that the photosphere is close to the NS surface and, thus, 
that the edge is gravitationally redshifted. For typical NS param- 
eters (e.g., M = 1.4M , R = 10 km) z ~ 0.3 at the surface and 
the intrinsic energy of the 8.5 keV edge would be 1 1 keV. This 
is close to both the hydrogen-like and helium-like edges of Ni 
(10.8 keV and 10.3 keV, respectively). That the edge may be due 
to Ni is supported by calculations which show that at the large 
ignition deptf0 implied by the intermedia te duration of the burs t 
from 4U 0614+09 (y base ~ 10 10 g crrr 2 ; iKuulkers et al] |2009). 
Ni is the dominant prod uct of constant pressure helium burning 
(lHashimotoetal.1119831) . 

The 4U 1722-30 moderate expansion phase edge at Sedge ~ 
1 1 keV is also close to the unshifted hydrogen-like Ni edge, 
consistent with the large measured photospheric radius during 
this phase. After the moderate expansion phase, when the pho- 
tosphere is at the NS surface, £ e dge ~ 7 ± 1 keV, consistent with 
a gravitationally redshifted hydrogen-like Ni edge. 

We can also use the measured optical depth of the edge, 
T opt , to demonstrate that not only must the absorption be due 
to a heavy element such as Ni, but that its abundance must be 
much higher than the solar value. We do so by assuming that 
T op t = Tbf{E e ) = N<Tbf(E e ), where N is the hydrogenic col- 
umn density of the element above the photosphere, <T\,f(E) st 
6.3 x 10~ 18 (£edge/£) 3 Z -2 cm 2 is the bound-free cross-section, 
and we take £ e dge = 13. 6Z 2 eV (see, e.g.. lBildsten et al] l2003). 
We approximate the hydrogenic column density as N(T,p,Z) = 
f(HN e , where / is the fraction of atoms belonging to that ele- 
ment (the 'number fraction'), is the fraction in the hydrogen- 
like state at a given temperature and density from Saha equilib- 
rium (we assume an element is either fully ionized or hydrogen- 
like), N e a* cr^} is the electron number column density, and <x T h 
is the Thomson cross-section. Such an approx imation should 
give a reasonable estimate of N (see, e.g., IWein berg et al. 
l2006bh . 

In the top panel of Fig. [5] we show as a function of effec- 
tive temperature kT s s and r p h = {L^&l AncrT^) 1 1 1 for 44 Ti and 
56 Ni at a densit)Q of p — 0.01 g crrr 3 . In the bottom panel we 

4 If instead it was a radial free-fall, the redshift could be significant 
even at these large radii. The ratio of observed to emitted frequency 
would be v obs /v em = 1 -/J where p = v/c = (2GM/c 2 r) 1/2 is the free fall 
velocity. Thus, for r - 30 km have v bs - 0.6v cm . 

5 We choose to express depth in terms of the mass contained in a 
column of unit cross-sectional area. 

6 In the wind models, dlnp/dhi r i -3 and p» lg cirr 3 at r = R 
(seejPaczvnski & Proszvnski 1986); since the observed expansion fac- 
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Fig. 5. Fraction of 44 Ti and 56 Ni in a hydrogenic state £h (top 
panel) and the optical depth of the edge T\,i(E e ) = Ncrbf(E e ) 
(bottom panel) as a function of effective temperature (bottom 
axis) and blackbody radius r p h = (L^dilAncrT^) 1 ^ 2 (top axis) 
assuming a density of p — 0.01 g cirT 3 . Results are shown for 
mass fractions of X — 0. 1 and 1 . 



show the implied optical depth of the edge T\,f(E e ) = No-u(E e ) 
for mass fractions of X = 0.1 and 1. As the figure shows, a mass 
fraction X 0. 1 is required in order to obtain an optical depth as 
large as the measured values T opt <: 1 at the observed tempera- 
tures. Furthermore, elements Z $ 22 cannot produce sufficiently 
large optical depths even for X = 1 . A Ni (Fe) mass fraction of 
X ~ 0.1 is ~ 1000 (100) times larger than the solar abundance. 

Rotational broadening is expected to be small compared to 
the PCA spectral resolution and the width of an absorption edge. 
The NS in 4U 0614+09 ha s a spin frequency of v spin = 415 
Hz (IStroh maver et al] 120071) . resulting in a maximum relative 
Doppler shift of v/c = 2n , v sp i n /?NS l c — 0.1 and a FWHM of the 
broadening profile approximately a factor 2 lower (n.b., viewing 
angle effects and the smaller rotation velocity of an expanded 
photosphere tend to further reduce the broadening). The other 
two bursters have no measured spin frequencies, but the high- 
est confirmed spin fr equency ever mea sured for a burster is only 
50% larger (620 Hz; iMuno et alJl2002h . 

6.2.2. Small radii during moderate expansion phase 

Steady state models of radiative winds from NSs predict photo- 
spheric radii of 100 - 1000 km for mass outflow rates of M » 
10" - 10 19 g s" 1 dPaczyriski & Pr6szvhskilll986t iJoss & Meli <3 
Il987t iNobili etafll 19941) . The measured photospheric radii dur- 



ing the moderate expansion are much smaller than this (r p h ~ 
30 - 50 km). The models were not calculated for lower M for nu- 
merical reasons and perhaps the observed radii are small because 
the outflow rates are M < 10 17 g s -1 . This seems unlikely, how- 
ever. From energy conservation. iPaczvhski & Pr6szvriskll (fl986) 
showed that the radiatively-driven mass outflow rate from a neu- 
tron star is given by M st aMo, where a = (Lbase - ^EddV^Edd is 
the excess luminosity at the base of the ignition layer relative to 
Eddington and Mq — 1.7 x 10 18 g s~' is a constant that depends 
only on the neutron star mass and radius (taken to be M = 1 .4M 
and R — 10 km). Such a low M (i.e., a 5 0.05) would therefore 
seem to require fine-tuning as it would imply that all the bursts 
are super-Eddington by only the slightest amount. 

An alternative possibility is that the outflow structure is trun- 
cated at a radius where absorption becomes sufficiently more 
important than electron scattering; material above this radius 
has a lower effective Eddington limit an d is blown away. The 
wind models either ignored a b sorptio n (Ebisuzaki et al. 1983; 
lPaczvriskT & Proszvriski Il986t INobili et al] 19941) or~ass umed 
the wind consists of solar abundances (IJoss & MeliaH l987). The 
presence of heavy element ashes in the wind can, however, sig- 
nificantly modify the structure of the outflow. 

To see why, note that at the temperatures observed during 
the moderate expansion phase (kT^ ^ 1 - 2 keV), the fraction 
of Ni in the hydrogen-like state is 0.001 - 0.4 from Saha equi- 
librium at a density p = 0.01 g crrT 3 (this assumes the Ni is 
either fully ionized or hydrogen-like; Ni is the dominant prod- 
uct of burning at these depths, see § 16.2. j} , The radiative force 
on hydrogenic Ni in the ground state is, in the optically thin 
limit (i.e., ignoring line blanketing) F Tal \ = nhe 2 f\^iE '(E) / 'm e c 2 
(BildstenetaL|2003; Mic haudl 19701) . where the flux per unit en- 
ergy at the Lya transition F(E) is nearly constant across the line 
(E = 8.0 keV) and fi_>7 = 0.4 2 is the oscillator strength for 
hydrogenic Ni (Fuh r et al.ll 1 98 lb . At the observed temperatures, 
the radiative acceleration on Ni is then a m d a 10 15 - 10 16 cm s~ 2 
assuming a blackbody spectrurrfl a factor of 100 larger than 
the gravitational acceleration at the observed photospheric radii 
r p h ~ 30 - 50 km. The force acting on the bound electrons is 
therefore ^ 100 times larger than that acting on the free elec- 
trons and line driving will become dynamically significant once 
the abun dance ratio of bound to free electrons exceeds ~ 1/100 
(see, e.g.. lGavlevll995l) . This suggests that once material reaches 
» 50 km and a significant fraction of ions become hydrogenic, 
line-driving will accelerate the flow to velocities much higher 
than the u « 10 3 km s found in the aforementioned wind mod- 
els. 

Since p a in a steady state wind, the structure of the 
outflow above « 50 km may be much more tenuous than the 
wind models suggest. If line-driving accelerates the material to 
~ 10 4 - 10 5 km s , the column density of overlying material 
will be a factor of ~ 10-100 smaller than that found in the 
wind models. The radius of the electron scattering photosphere 
may therefore be smaller by the same factor, in better agreement 
with the observed r p h. One caveat is that it is not clear whether 
lines or electron scattering will be the dominant source of opac- 
ity; a wind calculation that accounts for the absorption opacity 
from heavy elements is therefore needed in order to reliably es- 
timate the location of the photosphere. While the electron scat- 



tors are « 3 — 5 during the moderate expansion phase, we expect the 
photospheric density to be as 1CT 2 gem" 3 . For £ H £ 0.1, the optical 
depth Tbf scales almost linearly with p. 



1 The spectra of t ype I bursts are somewhat harder than a blackbody 
at the observed r, ff dLondon et al . 1986; Pa vlov et al .11 991; Mad ei et alj 
l2004l : lMaiczvna et al.ll2005r) . The radiative force acting on hydrogenic 
ions may therefore be larger than our estimates, although the ionization 
fraction will also be smaller. 
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tering opacity dominates the free-free and bound-free opacity 
for the measured temperatures and densities of the photosphere 
(even if the compo sition is primarily heavy elements; see, e.g., 
ICox & Giuh1 [T968). the bound-bound opacity may be important. 

6.3. Superexpansion duration 

If superexpansion is due to a super-Eddington flux driving the 
photosphere to large radii, as the time profiles suggest, why is 
the superexpansion duration f se always short (a few seconds) and 
independent of the decay time Td eC ay ~ 10 - 1000 s (see § [5])? 
Since the longer a burst, the longer the flux remains near its peak 
value, one might have expected f se to increase with increasing 
Tdecay, as is the case with the moderate expansion duration f me . 

We now describe a possible explanation for this lack of cor- 
relation between f se and Tdecay In Appendix ICl we show why the 
lack of correlation is unlikely the result of either a brief (second- 
long) phase of highly super-Eddington flux or hydrodynamic 
burning. 

X-ra y burst wind models have all assumed a steady state 



outflow (lEbisuzaki et al. 19831: Paczvriski & Proszvns 



idy state 



Uoss & Melialll987HNobili et al.ll 1994 and do not therefore de- 
scribe the time-dependent response of the overlying layers in the 
very first moments of radiative driving. These models do, how- 
ever, show that it takes « 1 s for tim e-independent con ditions to 
be established (see f c h ar in table 1 of lJoss & Meliain~987h . Given 
that f se » 1 s, this suggests that the superexpansion may be the 
result of a transient stage in the development of the wind. 

Analyzing this possibility requires a time-dependent wind 
calculation and is beyond the scope of this paper. It is nonethe- 
less worth noting that the superexpansion is consistent with the 
ejection of a geometrically thin shell of material of initial column 
depth i/ej s> f/ph, where j/ p h — 1 g crrT 2 is the depth of the photo- 
sphere on the NS surface. This shell, which may correspond to 
the layer of material above where the flux first becomes super- 
Eddington^, will become transparent when it reaches an expan- 
sion factor r/R ~ (j/ej/j/ph) 1 ^ 2 - Once it becomes transparent — at 
a time f se according to this picture — the observer suddenly sees 
the underlying neutron star again. Since the wind has had enough 
time to reach a steady state by this time, the photosphere is lo- 
cated at the moderate expansion value r p h ~ 30 - 50 km. 

The nuclear energy release from He burning 
(» 1 MeV nucleon ) is approximately 1% of the gravita- 
tional binding energy at the NS surface, and therefore at 
most 1% of the accreted mass can be ejected. This implies 
j/ej ~ 0.01i/base- If we assume the velocity of the ejected 
material is roughly constant with radius, the initial depth of 
the ejected column is y e j = y p h{v P htse/R) 2 and we have the 
constraint i/base ~ 10 8 (u p hf S e/10 4 km) 2 g cirT 2 . Since even short 
duration X-ray bursts have ignition depths i/base > 10 8 g cm -2 , 
this constraint is easily satisfied in the case of 4U 1722-30 
(u p h ~ 10 3 km s and t se — 5 s; see § 13.21 and Table [1). In the 
case of 4U 0614+09, the only other superexpansion burst for 
which there is an estimate of u p h, we get the strong constraint 

J/base ^ 10 9 g cirT 2 for i> p h ~ 3 x 10 4 kms _1 and f se = 1.2 s. 
This burst is an intermediate-duration burst and the ignition 



During the burst rise, the flux first exceeds the Eddington limit at 
some depth deep below the photosphere. This is because the electron 
scattering opacity increases towards th e surface due to th e temperature- 
dependent Klein-Nishina corrections (Paczyriski 1983). A given flux 
can therefore be sub-Eddington in the deeper layers and at a certain 
point become super-Eddington as it diffuses upward into lower density, 
cooler material. 



depth inferred from the cooling time is i ndeed consistent with 
i/base ~ 10 10 g cirr 2 dKuulkers et alJl2009h . 

7. Summary and conclusions 

The volume and quality of data on superexpansion bursts has in- 
creas ed substantially since th e first detailed study nearly 20 years 
ago (I van Paradiis et al.l[l99 0). We found 28 new superexpansion 
bursts in the burst catalogs and literature, several of which were 
detected with much higher temporal and spectral resolution than 
the 4 original superexpansion bursts. In studying this larger sam- 
ple, we found that superexpansion bursts have the following fea- 
tures: 

1 . At least 3 1 of the 32 superexpansion bursts are from (candi- 
date) ultracompact X-ray binaries. 

2. Three PCA bursts show significant spectral deviations from 
the absorbed black body model that typically describe X- 
ray burst spectra. The spectral deviations can be explained 
as absorption edges, with edge energies and depths that vary 
with time. 

3. Superexpansion is seen in short duration bursts, 
intermediate-duration bursts, and superbursts and the 
superexpansion phase is always followed by a moderate 
expansion phase. 

4. The duration of the superexpansion is always a few seconds, 
independent of burst duration, while the duration of the mod- 
erate expansion phase is proportional to burst duration. 

5. The photospheric radius during the moderate expansion 
phase is significantly smaller than that predicted by models 
of radiation-driven X-ray burst winds. 

We showed that the absorption edges and small photospheric 
radii of the moderate expansion phases may indicate the pres- 
ence of heavy-element ashes in the wind. In particular, the edge 
energies and optical depths appear consistent with a redshifted 
hydrogen-like or helium-like Ni edge and high Ni mass fractions 
(X ^ 0.1); Ni is expected to be the dominant product of burning 
at the large ignition depths inferred from the duration of these 
bursts. While we do not have a good explanation for why the 
edge energies seem to shift over the course of a few seconds, we 
speculate that it may be due to variations in which elements and 
ionization states dominate as r p h decreases and the temperature 
increases. 

The brevity of the superexpansion and its lack of correlation 
with burst duration may be the result of a transient stage in the 
wind's development, perhaps during which a shell of material is 
ejected to large radii by the sudden onset of super-Eddington flux 
deep below the photosphere. A time-dependent wind calculation 
using realistic compositions and opacities is needed in order to 
assess the viability of this proposed explanation. 

The spectral deviations seen in the PCA bursts suggest that 
the spectra of superexpansion bursts have the potential to pro- 
vide great insight into burst and NS physics. High spectral res- 
olution observations with sub-2 keV energy coverage, as can 
be achieved with Chandra or XMM-Newton, are clearly needed. 
The shortest burst recurrence times for X-ray sources that ex- 
hibit almost exclusively superexpansion bursts ( e.g., 4U 1722-30 
and 4 U 1812-12) are of order 5 to 10 days (e.g. Jin 'tZand et al.l 
2007), implying required exposure times of ~ 500 ks if one were 
to plan a dedicated observation. There is oft en some predictabil- 
ity to when an X-ray burst will occur (e.g JUbertini et alJll999t 
Gallowav et al. 2008a) and with careful planning less demand- 
ing exposure times (say by a factor of two) may be feasible. 
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Appendix A: PCA spectral data preparation and 
analysis method 

Details about the PCA data collecting modes, analysis settings, 
and time and spectral resolution are provided in Table IA.1I For 
4U 0614+09 and 4U 1722-30 data were obtained in the burst 
catcher mode. For 4U 1820-30 standard-2 data were obtained. 
The PCA spectra were accumulated from all active PCUs. While 
the spectra for 4U 1820-30 were drawn from only the upper 
xenon layer ('LR1'), those of the other two bursts were from 
all layers (the 'burst catcher' mode does not allow separation of 
PCUs or layers). 

A varying time resolution was employed in the time-resolved 
spectroscopy, based on the resolution available, the statistical 
quality of the data and the changes in the spectrum. The spectra 
were rebinned so as to have at least 15 photons per bin in or- 
der to ensure applicability of the^- 2 statistic. Contributions from 
particle and cosmic backgrounds, as predicted through version 
3.6 of the tool pcabackest, were subtracted. 

Response matrices were generated with version 11.7 of the 
tool pcarmf. A systematic error of 1 .5% was assumed per photon 
energy channel. The instrument team recently suggested^ that 
the systematic error should be set to 0.5%. While our fitted pa- 
rameter values do not change if we adopt their prescription, do- 
ing so increases the x 1 values somewhat. Perhaps the larger sys- 
tematic uncertainty is due to a spectrum (a few-keV black body) 
that is different from that of the Crab source (a power law with 
a photon index of 2. 1) which was solely used for the calibration; 
this should be investigated further. 

Burst spectra were corrected for detector dead time follow- 
ing the procedure described in the RXTE Cook BookFl. The time 
resolution for such a correction is limited to 0.125 s, so interpo- 
lations were necessary for the precursor in 4U 0614+09. 

The burst spectra were compared with mo dels between 3 
keV and 30 keV using XSPEC version 12.5.0 (lArnaudl fl996). 
The 3 keV threshold signifies the photon energy below which 
the instrument response is not well calibrated. The upper thresh- 
old signifies the photon energy above which the burst emission is 



9 http://www.universe.nasa.gov/xrays/programs/rxte/pca/ 
doc/rmf/pcarmf- 11.7/ 

10 http://heasarc.gsfc.nasa.gov/docs/xte/recipes/pca_deadtime.html 



in 't Zand and Weinberg: Superexpansion X-ray bursts 
Table A.l. Details of the data (analysis) for the 3 PCA superexpansion bursts. See also Figs. [2] [3] and lB.ll 
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4U 0614+09 


4U 1722-30 


4U 1820-30 1 


Start time in spacecraft clock seconds 2 


223941163.45 


90054034.88 


179459214.38 


UTC 


4-Feb-01 


8-Nov-96 


9-Sep-99 




21:52:34 


7:00:31 


1:47:58 


Active PCUs 


0,1,2,3 


0,1,2,3,4 


0,2,3 


Time interval of superexpansion 3 


0.04-1.2 


3.63-5.26 


15-17 


End time of moderate expansion 3 


89 


23 


1400 


Time interval with available burst catcher mode data 3 


0.93/35.93 


-0.25/3.25 


n/a 






7.75/97.50 




Time/channel resolution burst catcher mode data 


8ms/64 


2ms/64 


n/a 


Time/channel resolution science event mode data 


125 yus/64 


16yUS/64 


n/a 


Time/channel resolution resolution standard-2 data 


16 s/129 


16 s/129 


16 s/129 


Time interval of initial spectral study (Fig.IB.lt 3 


7.0-23.0 




325.1-1861.1 


Time interval of 'pre-burst' spectrum 3 


-740/- 100 


+5760/+9220 4 


-3140/-2320 



'4U 1820-30 data is accumulated only for the upper xenon layers; 2 In spacecraft clock seconds (i.e., seconds since Jan. 1, 1994, 0:00:00.00 UTC); 
3 Time in seconds since start time; 4 For 4U 1722-30 data were taken after the burst because pre-burst data are contaminated by a recent emergence 
of the spacecraft from the South Atlantic Geomagnetic Anomaly. 



expected to be negligible. At those energies, the spectrum is ex- 
pected to be dominated by the accretion flux. Thus, some handle 
is obtained on the accretion emission. 

Accounting for the changing accretion emission is not trivial. 
In super-Eddington bursts, one cannot simply subtract the spec- 
trum accumulated before the burst because the super-Eddington 
flux and expanding photosphere can affect the accretion disk 
(e.g., Ballantvne & Strohmaver 2004). Indeed, as shown in the 
main text, the observed flux drops below the pre-burst (i.e., ac- 
cretion) flux during the superexpansion phase. The flux during 
the moderate expansion phase may also affect the accretion. It 
is difficult to differentiate accretion flux from burst flux, because 
the burst emission dominates over the most sensitive part of the 
bandpass. Only at the highest energies does the accretion emis- 
sion dominate, but the statistical quality is insufficient to infer 
a spectral shape. We assume that the spectral shape of the ac- 
cretion emission remains, throughout the burst, identical to that 
of the pre-burst emission and fit its ('accretion') normalization. 
When possible we determined an empirical 3-30 keV spectral 
shape from data in a ~1000 s time interval prior to the burst (this 
is not always exactly 1000 s since there is not always data avail- 
able for a continuous 1000 s interval, see Table [AT). The data 
from 4U 0614+09 are consistent with a power law, and the data 
from 4U 1722-30 and 4U 1820-30 are consistent with a comp- 
tonized spectrum, all with low-energy absorption. This accretion 
model was incorporated in to our full spectral model of the data. 

Quoted errors on spectral parameters are single-parameter 
lcr values. 

Appendix B: PCA spectral model 

To understand the X-ray burst emission process, we first accu- 
mulated an optimum-quality spectrum for each burst by iden- 
tifying long data stretches when the spectral shape, as mea- 
sured by the flux ratio between two photon energy ranges, re- 
mains approximately constant. The times over which we accu- 
mulated data for this initial spectral study are given in Table lA. fl 
We then tested whether the resulting spectra are consistent with 
the canonical burst model, consisting of black body radiation 
(leaving free the temperature and normalization), low-energy 
absorption by cold gas of cosmic abundances (following the 
model for atomic cross sec tions and cosmic abundances of 
M orrison & M cCammonl ll983t leaving free the equivalent H- 



atom column density Nu) and accretion emission (leaving free 
the accretion normalization; see above). The fractional devia- 
tions of the spectra with respect to this model are shown in the 
upper row of Fig. lBTI Not shown are the results between 20 and 
30 keV because of the large error bars in that bandpass. None of 
the spectra fit the model. The shapes of the deviations are gen- 
erally similar with a maximum at 6 to 7 keV and a minimum 
between 10 and 13 keV. The amplitude of the deviations of the 
different bursts are similar and far larger than instrumental un- 
certainties (0.5% to 1.5% per channel). 

One approach to modeling the deviations is to include both a 
broad Gaussian emission line centered between 6 and 7 keV and 
an absorption edge close to 10 keV. A more physically moti- 
vated approach is to account for scattering ('reflection') and rel- 
ativistic blur ring of the burst's black body radiation by a hot ac - 
cretion di sk ( Ballan tvne & Strohmaverll2004t IBal lantvne 2004). 
iBallantvne & Strohmaverl d2004l) used such a reflection model to 
fit the deviations in the 4U 1820-30 superburst. Their reflection 
model accounts for reprocessing of thermal photons from the 
neutron star by the accretion disk (including comptonization, flu- 
orescence, recombination, and photo-electric absorption), which 
is allowed to come into thermal and ionization balance. The 
model includes a number of ionization stages of C, N, O, Mg, Si 
and Fe, and the parameters include the bolometric flux of the re- 
flection, the temperature of the black body, metal abundance and 
the ionization degree £ = 4nFx/nu of the scattering medium, 
with Fx the bolometric flux of the direc t black body rad iation 
and nu the hydrogen number density. See iBallantvnd (12004]) for 
further details of this model. 

We used the Ballantvne (2004) reflection model. In our im- 
plementation, we assume solar abundances and fit for two of 
the reflection parameters: £ and the bolometric flux. Assuming 
solar abundances is only relevant to Fe since it is the only el- 
ement in the model with spectral features in the PCA band- 
pass. The temperature is given by that of the black body com- 
ponent. Note that this model is quite different from reflection by 
a power law as commonly applied in AGN and Galactic black 
hole spectra. In bursts, the radiant energy per frequency decade 
(the 'vFy spectrum) peaks for black bodies at photon energies 
where photoelectric absorption by metals is very important, un- 
like for power law emission in AGN and Galactic black holes 
(Ballantvne 2004). 
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Fig.B.l. Spectral fits to the PCA bursts from 4U 0614+09 (left plots), 4U 1722-30 (middle plot) and 4U 1820-30 (right plots) 
during time intervals when the absorbed black body plus reflection model provides a good fit to the data (see Fig.|2]for example time 
intervals when the black body plus reflection model fails to provide a good fit to these bursts). The upper row shows the fractional 
deviations for a model that includes the absorbed black body component but not the reflection or absorption edge components. The 
lower row shows the deviations for a model that includes the absorbed black body and reflection component but not the absorption 
edge component. 4U 1722-30 is not shown because at that time does a reflection component explain the deviations. The spectra are 
taken during the moderate expansion. The error bars include the 1.5% systematic error. 



The reflection model also includes a convolution function for 
relativ istic blurring by the accretion disk following Fabia neTaT] 
( 1989). It is parametrized by the power law index for the radial 
dependence of the emissivity (set to -3), the inner disk radius 
(set to 10GM/c 2 with M the neutron star mass), the outer disk 
radius (set to 200GM/c 2 ) and the inclination angle (set to 30°). 
Whenever these parameters were left free during fitting, the im- 
provement in the fit was marginal. For practical purposes they 
were, therefore, kept fixed. 

The complete applied model can be written in XSPEC 
as: constant*(wabsi*(bbodyrad + rdblur*atable{reflection}) + 
wabs2*comptt/po) with constant the detector dead time factor, 
wabs\ and wabsi low-energy absorption (different for the two 
instances in this formula), bbodyrad the Planck function, rdblur 
the relativistic blurring function, reflection the model for the re- 
flecti on of the Planck function against the accretion disk follow- 
ing iBaflanT^n^OOO^coOTpff a comptonization component fol- 
lowing Titarchuk ( 1994) and po a power law function. The latter 
two components are representative of the non-burst emission. 
The model has six free parameters: Nh, k7bb, Rbb, the flux and 
£ of the reflection component and the accretion normalization of 
comptt ox po expressed relative to the pre-burst value. 

Including reflection in the model is highly successful in re- 
ducing the deviations in the spectra of 4U 0614+09 and 4U 
1820-30 for the particular time interval specified in Table IaTI 
(the "time interval of initial spectral study"). We find^- 2 =2.3 and 
0.75, respectively (the value for 4U 1820-30 is consistent with 



iBalla ntvne & Strohmaverl [2004l) . The remaining deviations are 
shown in the lower row of Fig. IB. II The black body plus reflec- 
tion model does not, however, fit all the spectra of 4U 0614+09 
and 4U 1820-30 (see § 13. It . Furthermore, it almost never pro- 
vides a good fit to the spectra of 4U 1722-30 (the middle-right 
panel in Fig.|2]is representative of the typical spectral deviations 
found at all times in black body plus reflection model fits to the 
4U 1722-30 spectra). 

Even when the black body plus reflection model provides a 
good fit, the flux of the reflection component, a free parameter, is 
often higher than that of the direct black body component. This 
is due to the fact that the fitting procedure primarily tries to fit 
the narrow features - the broad emission line and the absorption 
edge. Since the continua of the reflection and direct black body 
components are very similar in the PCA bandpass, their contri- 
butions are di fficult to disentangle. High reflec tion fluxes were 
also found by Balla ntvne & Strohmaverl (120041) in the late-time 
burst spectra of 4U 1820-30 (see their Fig. 1). As a result of a 
high reflection flux, the flux of the direct black body component 
becomes so low that the black body radius becomes substan- 
tially lower than 10 km. This is physically unrealistic. Without 
a proper model for the geometry of the accretion disk and inde- 
pendent measurements of the iron abundance in the donor star, it 
is impossible to meaningfully constrain the reflection flux in an 
independent manner. 

We calculated bolometric luminosities and emission region 
radii by fitting the spectra without the reflection (or edge) com- 
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ponent. After fitting, the absorption and the accretion normal- 
ization was set to zero, the deadtime correction set to 1, and the 
response matrix in XSPEC set to a dummy value. The flux could 
then be estimated in a broader energy range than the PCA band- 
pass. The unabsorbed flux was calculated between 0.05 and 30 
keV, essentially a bolometric bandpass for a black body with a 
temperature between 0.5 and 3 keV. Since the fits are often for- 
mally unacceptable, no errors were calculated. 

We tried several other functions to model the spectra. These 
include combinations of two different black bodies, combina- 
tions of one black body with a power law (with a free spec- 
tral index) , gauss function or relativistic iron line (following 
lLaorlll99lL model laor in XSPEC), and single reflection mod- 
els of a power law against a neutral or ionized disk (following 
Magdziarz & Zdziarski 1995, models pexrav/pexriv in XSPEC). 
While some of these models fit some of the spectra reasonably 
well, they only do so over limited time intervals in any given 
burst. A broken power law fits the data with the same goodness 
(xl) at any time, but there is no physical basis for such a model. 



Appendix C: Ruling out alternative explanations for 
the lack of correlation between t se and T deC ay 

Here we first consider the possibility that the transition between 
the superexpansion phase and the moderate expansion phase is 
due not to a shell ejection, as suggested in § 16.31 but rather to 
a transition in the luminosity at the NS surface. In particular, if 
the luminosity exceeds not only the local Eddington luminosity, 
L Edd = AjiGMc/kq, but also the Eddington luminosity at infinity, 
^Edd.oo — iEdd(l + z) 2 , a strong radiation-driven wind will form 
(/fo = 0.2 cm g is the hydrogen-free electron scattering opac- 
ity). Such a wind can potentially drive the photosphere to very 
large radii (k, 100 km). If, however, Z^dd < L < LEdd.oo, a com- 
paratively mild quasi-static expansion is su fficient to reduce the 
luminosity at the photo s phere below L E dd (Hanawa & Sugimoto 



20 



luminosity at the photosphere below L E dd (Hanawa & bugimoto 
Il982t lPaczvhskilfl983l: iPaczvnski &"P r6szvnski Tl986t) . In this 
case, the PRE will be less extreme. This suggests that perhaps 
the large radius expansion factors reached during the superex- 
pansion occur because initially L > Z,Edd,°o- Then, after a time 
f se « seconds, the luminosity decreases to Z^dd < L < Z,Edd,c» 
and the photosphere retreats closer to the NS surface, marking 
the start of the moderate expansion phase. This phase continues 
for a duration of order the cooling time after which the luminos- 
ity decreases to L < Z^dd and the photosphere settles on to the 
NS surface. 

To test this idea, we carry out time-dependent cooling calcu- 
lations for a range of ignition depths. To s olve the heat transport 
equati o n we use the technique describ ed in lWeinberg & Bildstenl 
(2007j); ICumming& Macbeth] d2004 carry out similar calcula- 
tions for superbursts. We assume that the atmosphere is radia- 
tive and initially at its peak temper ature, in accord with models 
of the rise of X-ray bur sts (see e.g.,|Cumming & Bildsten 2000; 
IWeinberg et afl l2006b). We assume a peak temperature at the 
base of the ignition layer i/ base of r base ,initiai = /T max , where / is 
factor less than but of order unity and T m - dx = Qgi/base/a) 1 ^ 4 is the 
critical temperature at which radiation pressure completely dom- 
inates (g = ( 1 + z)GM/R 2 is the gravitational acceleration and a 
is the radiation constant). We chose the value of / by solving for 
Tbase.initiai such that the fluence (the cooling curve integrated over 
the duration of the burst) equals the total nuclear energy release 
E = 47r/? 2 (/base2nuc/(l + z) assuming gnuc = 1.6 MeV nucleon , 
corresponding to helium burning completely to nickel. Although 
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Fig. C.l. Model calculations of the radiative cooling following 
the burst rise. The solid lines show the luminosity L for ignition 
depths i/ base = 3 x 10 s , 10 9 , and 10 10 g cirr 2 with / = 0.88, 0.85, 
and 0.78, respectively. The values for / were chosen such that 
the fluence equals the nuclear energy released in burning the en- 
tire He column to Ni. The Eddington limit at the NS surface and 
at infinity (assuming z = 0.31 and M — 1.4 M ) are drawn as 
dotted lines. 



our conclusions are not sensitive to the exact value of /, for ref- 
erence, we find / » 0.8 - 0.9. 

In Fig. IC.ll we show time profiles of the calculated luminos- 
ity L for three ignition depths i/base- The dotted lines mark Z^dd 
and Z,Edd,co for a 1.4 M NS and an opacity kq — 0.2 cm 2 g _1 . 
An observer would see an Eddington-limited light curve where 
L > L E dd (the excess luminosity is used to drive the outward ex- 
pansion of the overlying layers) and a light curve that follows 
the plotted luminosity where L < LEdd- 

As the figure shows, the profiles all have the same shape re- 
gardless of i/base and /: an early flat phase followed by a break 
and a steep decay. The factor / primarily determines the overall 
vertical scale while j/ base determines the thermal time of the igni- 
tion layer and therefore the location of the break. The key result 
is that the early-time light curve is very flat, barely decreasing 
until after the break, at which point it rapidly falls below LEdd- 
Such a light curve cannot explain the features of the superexpan- 
sion as there is no second-long, depth-independent phase during 
which L > Z,Edd,oo- 

The light curve is flat because we assumed that the initial 
profile was the radiative profile (T oc y n with n as 1 /4) and there- 
fore L oc T 4 /y is constant until the thermal wave reaches i/base 
and the ignition layer begins to cool. If we had assumed a shal- 
lower initial profile, the early-time flux w ould decrease with time 
(see, e.g., the superburst calculations of Cumming & Macbeth 
2004 who assume an initial profile with n — 1/8 corresponding 
to an isobaric deflagration). However, absent any strong motiva- 
tion for assuming a non-radiative initial profile, we conclude that 
a variation in the early time L is unlikely to explain the lack of 
correlation between f se and Tdecay 
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in 't Zand and Weinberg: Superexpansion X-ray bursts 



Since the superexpansion phase always occurs within sec- 
onds of burst onset, perhaps it is associated not with post-peak 
cooling, as assumed above, but rather with the physics of the 
burning during the burst rise. Due to their larger ignition depths 
and helium mass fractions, the peak energy generation rate in 
superexpansion bursts can be considerably higher than in ordi- 
nary short-duration bursts. If the ignition is deep enough, the 
heating time due to burning fheat = (dlnTb/df) _1 can become 
shorter than the local dynamical time fdyn ~ 10~ 6 s. Such hy- 
drodynamic burning could drive strong shocks through the up- 



2006aj). Ho wever, when we s olve the rise using the prescription 
described in I Weinberg et al.l (l2006bl) assuming an initially pure 
He layer, we find that f nea t > fdyn unless y\, ase <: 10 11 gem -2 . 
This suggests that strong shocks do not form during the rise of 
superexpansion bursts as they ignite at j/base ~ 10 10 g cm" 2 . 




